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Abstract 

We study that a minimal supersymmetric standard model with an extra U(l)' 
gauge symmetry may accommodate the explicit CP violation at the one-loop level 
through radiative corrections. This model is CP conserving at the tree level and 
cannot realize the spontaneous CP violation for a wide parameter space at the one- 
loop level. In explicit CP violation scenario, we calculate the Higgs boson masses 
and the magnitude of the scalar-pseudoscalar mixings in this model at the one-loop 
level by taking into account the contributions of top quarks, bottom quarks, exotic 
quarks, and their superpartners. In particular, we investigate how the exotic quarks 
and squarks would affect the scalar-pseudoscalar mixings. It is observed that the 
size of the mixing between the heaviest scalar and pseudoscalar Higgs bosons is 
changed up to 20 % by a complex phase originated from the exotic quark sector of 
this model. 
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I. INTRODUCTION 



Possessing two Higgs doublets is one of the characteristics of the minimal supersymmetric 
standard model (MSSM) in order to give masses to the up-quark sector and the down- 
quark sector separately [1]. This property enables in principle CP violation to occur in 
the Higgs sector of the MSSM through the mixing between the CP even and the CP odd 
states [2]. The mechanism of CP violation in the MSSM has been studied by many authors 
[3-12]. It is found that at the tree level, neither explicit nor spontaneous CP violation is 
possible in the Higgs sector of the MSSM, because any complex phases in it can always 
be eliminated by rotating the Higgs fields. Even at the one-loop level, spontaneous CP 
violation is disfavored because it requires a very light neutral Higgs boson, which has 
already been ruled out by experiments [3,4]. On the other hand, explicit CP violation is 
viable in the MSSM at the one-loop level by virtue of the radiative corrections due to the 
loops of relevant particles, such as quarks and squarks [5-12]. The radiative corrections 
by these particles yield the mixing between the CP even and the CP odd neutral Higgs 
bosons. Thus, it is possible to achieve the explicit CP violation in the radiatively corrected 
Higgs sector of the MSSM. 

However, a drawback of the MSSM is that it contains a Higgs mixing term HH1H2, 
where \i is a parameter having mass dimension, and Hi (i = 1, 2) are Higgs doublets. It 
is known that this [i parameter causes a hierarchy problem with respect to the symmetry 
breaking scale [13]. A natural framework to solve this problem is to generate the \x 
parameter via the vacuum expectation value ( VEV) of a Higgs singlet [14] . A possibility 
for the framework is to extend simply the Higgs sector of the MSSM by introducing a Higgs 
singlet, like the next-to-minimal supersymmetric standard model, where the \x parameter 
is practically replaced by a dimensionless parameter times the VEV of the Higgs singlet 
[14]. 

A more plausible possibility to solve the \x problem is to extend not only the Higgs 
sector but also the gauge symmetry of the MSSM by introducing a nonanomalous U(l)' 
symmetry broken at the TeV scale [15]. Extending the gauge symmetry to include an 
additional 17(1)' symmetry is widely considered in various theoretical models, such as 
string models or GUT models. The MSSM with an extra U(l)' symmetry can actually 
forbid the \1H\H1 term but allow \SHiH 2 , by assigning appropriate U(l)' charges, where 
S is the Higgs singlet and A is a dimensionless parameter. After S develops a VEV of the 
electroweak scale, this A term would effectively generate the \i term at the electroweak 
scale [16]. 

In the literature, a number of studies have been performed on the MSSM with an extra 
£7(1)' symmetry [17-20]. Several authors have calculated the Higgs boson masses from 
the tree-level Higgs potential of this model [17] as well as from the radiatively corrected 
Higgs potential [18,19]. It is found that, at the tree level, the lightest scalar Higgs boson 
in this model can be heavier than Z boson. That is, the upper bound on the mass of the 
lightest scalar Higgs boson in this model is larger than the one in the MSSM without an 
extra U(l)' symmetry. Moreover, it is also found that this tree-level mass of the lightest 
scalar Higgs boson in this model can be significantly affected by radiative corrections. 

Recently, this model has been further investigated within the context of the supersym- 



2 



metric CP problem [20], as it can have explicit CP violating phases in the Higgs sector. At 
the tree level, like the MSSM, this model has no CP mixing between the CP even and the 
CP odd states in the Higgs sector because the single CP phase arising from the tree-level 
Higgs potential can always be eliminated by rotating the Higgs fields. In Ref. [20], the 
one-loop contribution from the top and stop quark loops are considered in the explicit CP 
violation scenario. It is observed that the lightest neutral Higgs boson remains essentially 
CP even for a wide parameter space, with a fixed value of tan (5 — 1, where tan (3 = Vifv\ 
is the ratio of the VEVs of two Higgs doublets. Meanwhile, large CP mixtures for the 
other heavier Higgs bosons may be realized as the size of the effective /x terms becomes 
large. 

In this paper, motivated by the results of Ref. [20] , we would like to study this model 
in some detail in the explicit CP violation scenario. The radiative CP mixing may be 
generated by means of the complex phases coming from the one-loop effective potential 
due to top quarks, bottom quarks, the exotic quarks, and their superparticles, where the 
exotic quarks and squarks with electric charges ±1/3 are introduced into this model in 
order to cancel gauge anomaly [21]. In particular, we are interested in the contributions 
of exotic quark sector in this model. At the one-loop level, the exotic quarks and squarks 
are found to play some recognizable role in the scalar-pseudoscalar mixings. We find 
that the relevant CP phase from the exotic quark sector may change the size of the CP 
mixing between the heaviest scalar and the pseudoscalar Higgs bosons by up to 20 %. 
Also we find that this model yields a negative or extremely small axion mass, as well 
as other unacceptable predictions, at the one-loop level for full parameter space in the 
spontaneous CP violation scenario. Thus, we note that the present model is impossible 
to realize spontaneous CP violation at the one-loop level. 



II. HIGGS SECTOR 

GUT, such as the string-inspired E 6 model, might be a natural motivation to enlarge the 
gauge symmetry in order to accommodate an extra U(l)' symmetry into the MSSM. The 
E 6 gauge group may be decomposed into 

E 6 D SU(10) x U{1)^ D SU(5) x U(l) x x U(l)^ , (1) 

where £77(5) is further broken down to the Standard Model (SM) gauge group, SU(3)c x 
SU{2) L x U(l)y- At the electroweak scale, the desired extra £7(1)' symmetry may be 
given as an orthogonal linear combination of U(l) x and U(l)^ as 

Q' = cos6 E Q x + sin E Q , (2) 

where Q', Q x , and are the U(l)', U(l) x , and U(l)y charges, respectively. Sometimes, 
for four different values of the angle 9 E , the four particular combinations of the U(l) x 
and U(l)ip are called as follows: the x-model for 9 E = 0, the -^-model for 6 E = ir/2, the 

77-model for 9 E = tan _1 (— ^5/3), and the z/-model for 6 E = tan" 1 \/l5. The gauge group 
we consider is thus G = SU(3) C x SU(2) L x U(l)y x U(l)'. 
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The Higgs sector of the MSSM with an extra U(l)' symmetry consists of two Higgs 
doublets Hi = (i7°,i7f) and H 2 = (H^H®), and a neutral Higgs singlet S. This model 
possesses an extra pair of 577(2) singlet quarks, Dl and Dr, with electric charges —1/3 
and +1/3, respectively, introduced by the requirement of the gauge anomaly cancellation. 
For the fermion matter fields in this model, we take into account only the third generation 
of quarks, besides the exotic quarks. Then, the superpotential for the model we consider 
may be expressed as 

W « \SHiH 2 + h t QH 2 t c R + h b QHib c R + kSD L D R , (3) 

where Q is the left-handed quark doublet, t c R and b c R are the charge conjugate of the right- 
handed top quark and bottom quark, respectively, with h t and hb being their respective 
Yukawa coupling coefficients, and HiH 2 = H^H® — HiH 2 . All the coupling coefficients 
are dimensionless. This superpotential is symmetric under SU(3)c x SU(2)l X C/(l)y X 
U(l)', and the U(l)' symmetry is broken down by the VEV of S. 

The Higgs potential of this model at the tree level may expressed as a sum of the 
F-term, the D-term, and the soft breaking term, that is, 

V = V F + V D + V s , (4) 

where 

V F = \X\ 2 [(\Hi\ 2 + \H 2 \ 2 )\S\ 2 + \HiH 2 ^ 



'2| 

V D = 9 4{H\aHi + HlaH 2 ) 2 + 9 ^{\Hi\ 2 -\H 2 



2\2 



'■-> 



^(QilHil 2 + Q 2 \H 2 \ 2 + Q 3 \S\ 2 ) 2 

,2| tt |2 i ™2| tt |2 i ™2| c|2 



7 S = m 2 \Hi\ 2 + m 2 2 \H 2 \ 2 + m 2 3 \S\ 2 -[\A x (HiH 2 )S + H.c.} , (5) 

where (i = 1,2,3) are soft masses, g 2 , gi, and g[ are respectively the coupling coefficients 
of S77(2), U(l), and f/(l)'; B = (ai, a 2 , a^) are the Pauli matrices; Qi, Q 2 , and Q3 are 
respectively the U(l)' charges of Hi, H 2 , and S, satisfying the condition of Q1+Q2+Q3 = 
to obey the gauge invariance of the superpotential under £7(1)'. 

At the electroweak scale, the three neutral Higgs fields develop VEVs as < Hf >— Vi, 
< H® >= v 2 , and < S >= se 1 ^ 3 , where <p s is the relative phase between HiH 2 and S. In 
general, \A\ in the soft breaking term may be complex, like XAxe^ . However, it is always 
possible to make the phase <fi and the phase <fi s cancel out each other by redefining the 
Higgs singlet S, and thus the tree-level Higgs potential can be made free of any complex 
phase. Therefore, the CP symmetry is not violated in this model at the tree level. 

Now let us consider the one-loop radiative corrections to the tree-level Higgs potential. 
In supersymmetric models, the incomplete cancellation between ordinary particles and 
their superpartners yield the one-loop corrections to the tree-level Higgs boson masses. 
Generally, the most dominant part of the one-loop corrections to the tree-level Higgs 
potential come primarily from the top and stop quark loops. By considering the top 
quark sector only, it has been observed that explicit CP violation in this model is viable 
through the radiatively corrected Higgs potential [20]. For large tan/3, the contribution 
of the bottom and sbottom quark loops can also be significantly large. We would like to 
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consider in this paper the one-loop contributions from both the top and bottom quark 
sector, as well as from the exotic quark sector. 

To start with, we calculate the tree-level masses of the relevant particles in order 
to evaluate the radiative corrections to the neutral Higgs sector. The fermion matter 
fields obtain their masses after the electroweak symmetry breaking as m 2 = (h t V2) 2 for 
top quark, m\ = (h b Vi) 2 for bottom quark, and m\ = (ks) 2 for the exotic quark. The 
tree-level masses of their superpartners are given as 



m\ x i 2 = m\ + m\ =f m t \j A 2 + X 2 s 2 cot 2 (3 — 2XA t s cot (3 cos fa , 

m| g 2 = m\ + m\ =f m b \j A\ + X 2 s 2 tan 2 (3 — 2XA b s tan (3 cos fa , 

m !i k 2 = m K + m l T m k\j A\ + X 2 v 4 sin 2 2(3 /(As 2 ) - XA k v 2 sin 2(3 cos fa/s , (6) 

where tan/3 = v 2 /vi, v 2 = v \ + v\. Further, m T , m B , and m K are the soft SUSY 
breaking masses respectively for stop quarks, sbottom quarks, and the exotic squarks, and 
likewise A t , A b , and A k are respectively their trilinear soft SUSY breaking parameters of 
the mass dimension. One can notice that there are three phases <f) t , fa, and 0^ in the 
above expressions. These complex phases are determined by the generally complex A q 
(q = t, b, k) and the phase of s, the VEV of the Higgs singlet. We note that the D-terms 
do not contribute to the squark masses. 

The full Higgs potential at the one-loop level may be written as 

V = V + Vi, 

where V\ is the radiative corrections due to the relevant particles and their superpartners. 
According to the effective potential method [22], V\ is given as 

Ml 3- 



Vi = E n,Mt 



64vr 2 



log 



A 2 2 



where A is the renormalization scale in the modified minimal subtraction scheme, and the 
subscript I stands for various participating particles: t, b, k, t\, ti-, &i, 62, k\ and k 2 . The 
degrees of freedom for these particles including the sign convention are nt = n b = n k = 
— 12 and n~ t . = n~ h . = n~ k . =6 [i = 1,2), since in the above formula enter fermions with a 
negative sign while bosons with a positive sign. 

If CP violation takes place in the Higgs sector, the neutral Higgs bosons would not 
have definite CP parity, hence mixings among them. In case of explicit CP violation, the 
non-trivial tadpole minimum condition with respect to the pseudoscalar component of 
the Higgs field is given as 

3m 2 A t sin fa 2 2 ?>m 2 b A b smfa 2 2 

y 167r 2 v 2 sin 2 /3' /v h ' t2 ' \Qix 2 v 2 cos 2 [3 J V bl b2 ' 

3m 2 k A k smfa 2 2 
ttt, / {mi , m~ , 8 

\Q-jl2g2 J v fei' k 2 ' ' v ' 

where the dimensionless function / arising from radiative corrections is defined as 

1 



f(™ x > m y 



{m 2 y - m 2 ) 



2, m 2 x 2 m 2 
m x lo S -J2 ~ m y lo § 
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and the phase O is given by 4>o = 4> s + 0, with being the phase of \A\ in V and 4> s 
being the phase of s. Note that </>o is generally nonzero at the one-loop level whereas at 
the tree level it can be made zero by rotating the Higgs singlet. In the above tadpole 
minimum condition, the first term comes from the tree-level Higgs potential, and the 
remaining three terms come respectively from the top squark, bottom squark, and exotic 
squark contributions. 

At the tree level, the Higgs sector in the MSSM with an extra U(l)' has ten real 
degrees of freedom, which are decomposed by two neutral Goldstone bosons, a pair of 
charged Goldstone bosons, four neutral Higgs bosons and a pair of charged Higgs bosons. 
After the electroweak symmetry breaking, the two neutral Goldstone bosons and a pair of 
charged Goldstone bosons will eventually be absorbed into the longitudinal component of 
Z, Z', and W gauge bosons. Since the CP symmetry is conserved in the Higgs sector at 
the tree level, the four neutral Higgs bosons can be divided into three scalar Higgs bosons 
and one pseudoscalar Higgs boson, according to the CP parity. 

The squared mass matrix M of the four neutral Higgs bosons is given as a symmetric 
4x4 matrix that is obtained by the second derivatives of the Higgs potential with respect 
to the four Higgs fields, in the (hi, h 2 , h 3 , h&) basis. At the tree level, M is given by V° as 



~ Ml Ml Ml ' ® 
\ m\ ) 

where the tree-level squared mass of the pseudoscalar Higgs boson is given as 

2 QU c °s 00 
m A = 2XA x v , 

sin 2a 

with 

v 

tan a = — sin 23 

2s H 

standing for the splitting between an extra U(l)' symmetry breaking scale and the elec- 
troweak scale. Note that M° is divided into two blocks of the upper-left 3x3 submatrix 
and the single element M° 4 = m\, corresponding to the scalar part and the pseudoscalar 
part. There is no scalar-psuedoscalar mixing in M°, since Ml (i = 1 — 3) are zero, and 
hence no CP violation. Explicit calculations yield 

M ii = m z c °s 2 P + Zg'iQlv 2 cos 2 3 + m\ sin 2 3 cos 2 a , 

Ml = m% sin 2 3 + 2g[ 2 Q 2 2 v 2 sin 2 3 + ra\ cos 2 3 cos 2 a , 

Ml = 2g' 1 2 Q 2 3 s 2 + A 2 sm 2 a , 

Ml = g 2 QiQf} 2 sin 23 + ( A V - m| /2) sin 23 - m 2 A cos 3 sin 3 cos 2 a , 

Ml = 2g 2 QiQ 3 vs cos 3 + 2\ 2 vs cos 3 — m 2 A sin 3 cos a sin a , 

Ml = 2g 2 Q 2 Q3Vs sin/3 + 2X 2 vs sin/3 — m 2 A cos 3 cos a sin a , (10) 

At the one-loop level, the squared mass matrix may be decomposed as 

M = M° + M 1 
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where M° has exactly the same appearance as M° but m\ in M° should be replaced by 
m\ in M°, where m\ is given as 

2\v 



m , 



sin 2a 



377l|A t C08^t ,/ 2 

A A cos0 o - ^_ 9 „ 9 o f( m u 



3ml At, cos ^fe 
16tt 2 v 2 cos 2 /3 



167r 2 f 2 sin /3' 



6 2 > 



3m 2 .A k cos * 

167T 2 S 2 



! if 2 



rn 



fc 2 > 



11] 



One may notice that in the expression for m\, the first term comes from the tree-level 
Higgs potential and the remanning three terms comes from the contributions of the 
squarks. Therefore, M° contains not only the tree-level results but also the contribu- 
tions from radiative corrections. 

Now, let us calculate M 1 , which is obtained from V\. We may conveniently decompose 
M 1 as 

M 1 = M l + M b + M k 

where symmetric matrices M t , M b , and M k are obtained respectively from the top quark 
sector, the bottom quark sector, and the exotic quark sector contributions to V±, after 
imposing the tadpole minimum condition. 

Somewhat lengthy calculations yield the matrix elements of M t as 



M\ 2 



M, 



33 



M; 



14 



Mi 



12 



Mis 



M\s 



ML 



3m 4 A 2 s 2 A 2 g\ 



8tt 2 v 2 sin (3 (rn 2 2 — m 2 ^) 2 
3m\A\A\ g(m\, m\ t 
87r 2 ^ 2 sin 2 (3 {m 



3m; 



^ 8ir 2 v 2 sin 2 f3 
3m 4 A 2 A| g(ml, 

8tt 2 tan 2 j3 (m 

3m A t \ 2 A 2 s 2 sin 




3m A t A t A h log(m|/m|) 
4 7r 2 t; 2 s j n 2 p {m 2 ^ — m 2 ^) 



t 9^ m h m X> 



8tc 2 v 2 sin 4 P cos 2 a (m~ 2 — m 2 ^) 2 

_ 3m A t \A t sA h A h g{m\, m|) 

8ti 2 v 2 sin 2 P (m 2 ^ — m^) 2 
2 \ 



3m 4 A 2 sA? 



g{m\ m\) 3m 2 X 2 scosp 



3mt\sA ii log(m|/m|; 
8ti 2 v 2 sin 2 P (m 2 — m~ ) 



8tt 2 v sin/? tan /5 (m 2 ^ — m 2 ^) 2 87r 2 wsin P 



fH 



3m 4 A 2 A t s 2 A h sin t ^(m?, 



2 



t 2 > 



87r 2 f 2 sin P cos a 



(m| -m 2 h ) 2 
2 ^,2 \ 



8rc 2 v sin /5 tan P (m~ 2 



3m 4 Acos/3A fi log(m|/m|) 



m 2 r y 

tl ' 



3m\\A 2 t sA~ t2 sin 



it ^(H' m |) 



^7r 2 t> sin 2 P (m 2 2 



87r 2 v 2 sin P cos a 



+ 



3m 4 AAssin0 t log(m|/m|) 



m 2 — m 2 ) 2 ' 87r 2 w 2 sin 3 /3 cos a (m~ — m 2 ) 



t 2 
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3mfX 2 A t sA k sm^ t g(m\, m |) 
87r 2 f sin 2 f3 tan (3 cos a (m~ — m t 



JW 34 - o_9„.„;„2 «4.„„ «„„„ „ .1^1 „„2 \2 > 



where 



A^ = A t cos 0t — As cot (3 , 

A t ~ 2 = At — As cot j3 cos <f>t , (13) 
and the dimensionless function g is defined as 

2 2 , m 2 + m\ rn 2 
g(m x , m) = — § log — f + 2 . 

y — m x 

Likewise, the matrix elements of M b as 

3miA 2 A 2 bi g(m 2 bi , m|) 3m 4 A A^ log^/m^) 



M 



m: 



Ml 



ML 



M 



Mi 



M: 



where 
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s 7r 2 f 2 cos 2 /3 (m 2 — m 2 ) 2 Air 2 v 2 cos 2 /3 (m~ — m~ 

4 / to, 2 m. 2 



I 3m » log h b2 
8tt 2 v 2 cos 2 /? I m\ 

3mi\ 2 s 2 A 2 b2 g(m 2 bi , m|) 



87r 2 u 2 cos 2 /9 (m 2 — m\ 

y o 2 6 

3m 4 A 2 A? 3(m. 2 , m? ) 



33 «7T 2 cot 2 5 (m? -m? ) 2 

v o 2 Oi' 



3m 4 A 2 ^s 2 sin 2 b g(mf • 



fi' 



87r 2 f 2 cos 4 (3 cos 2 a (m~ — m~ ) 2 



3m 4 AA b sA 5i A 02 s(m|, mfj 3m 4 AsA, 2 log^/m 2 ; 



87T 2 i; 2 cos 2 /3 (m~ — m? ) 2 87r 2 i; 2 cos 2 (3 (m~ — m 2 b \ 
3m 4 AA 6 Ar At, , mfj 3m 4 A sin/3Ar log (m? /m 



87r 2 f cos /5 cot /5 (m~ — m~ ) 2 87r 2, ucos 2 /3 (m~ — m| ) ' 
, 3m 4 AA 2 sA^sin0 o ^(m^, m\) 3m\\A b s sin O log^/mfj 



87r 2 f 2 cos 3 /3 cos a (m~ — m 2 ) 2 8ir 2 v 2 cos 3 /? cos a (m~ — m 2 b \ 



3m\\ 2 sAl gjmlmj) _ 3m 2 A 2 stan/3 /(m2 ^ ^ 
7r 2 w cos /3 cot /5 (m~ — m~ ) 2 87r 2 t> cos /5 61 ' 62 ' 



6 3m 4 A 2 A fe s 2 A j)2 sin0 o g(m^, m|) 

87r 2 w 2 cos 3 3 cos a (m- — m~ ) 2 
3m 4 b X 2 A b sA- b2 sm(f) b g(m 2 h , m? 1 



M b = - - ^ 1 -2— , (14) 

34 87r 2 t> cos 2 /3cot /9 cosa (m? 2 — m^) 2 



Ag i = A b — As tan j3 cos </> , 

At = A b cos O — As tan j3 , (15) 



8 



and the matrix elements of M k as 



3ml\ 2 v 2 sin 2 /?A~ <7( m f > TO | 2 ) 

v fe2 fel/ 



M fc 



M: 



Mi 



M, 



M; 



3mf,\ 2 v 2 cos 2 /3A| 9(^1^ m\ , 



22 8tt 2 s 2 (m? - m? ) 2 ' 

v K2 fel ' 

33 87r 2 s 2 (m| - m| ) 2 47r 2 s 2 (m| - mf ) 

3m 4 . , ( ml ml 



8^5^ \ m k 



3m 4 fc A 2 ^ 2 sin 2 fc <?(<, mf^ 



8 7r 2 s 2 cos 2 a (m~ — ml ) 2 



3^A 2 , 2 sin2^Al i gfrn^ mfj 3 m 2 AV sin 2/3 
12 16tt 2 s 2 (m?-m?) 2 16vr 2 s 2 A *i' 



3mi\A k v sin (3A~ kl A~ k2 g{m\, m\ t 



13 8ttV (m 2 k2 - m 2 ki ) 2 

3m^sm/3A^ log^/m 2 ^) 
87r 2 s 2 (ml — m~ ) 

fel 

Mfc _ 3mtA 2 v4^ 2 sin/jA^sin0 fc 3(7^, 



87r 2 s 2 cosa (m| — ; 

3mj\A k v C osPA~ k A~ k2 g[m\, 

8ttV (mf -ml) 2 

K2 fel 

3m 4 A?; cos /3A^ log^/mfj 
87r 2 s 2 (m| — m| ) ' 
3mtA 2 A^ 2 cos/3Afe iS in0 fe #(mf i5 mfj 



24 87r 2 s 2 cosa (m\ — m\ ) 2 ' 

fe2 fel 

Mfc = 3m 4 AA 2 t;A^sin0 fc ^K i , m| 2 ) ^ 3mjXA k v sin fc logK^/mfJ 
87r 2 s 2 cosa (m~ — mf ) 2 87r 2 s 2 cosa (m~ — m~ ) ' 

fe2 fel fe2 fel 

where 

A£ = A k cos 0fc — Xv tan a , 

A^, 2 = A fe — Ai> tan a cos 0& . (17) 

Notice that the matrix elements M* 4 , Mf 4 , and M i4 (z = 1 — 3) are proportional to 
sin 4>t, sin 06, and sin0fc, respectively. If the three CP phases are not all zero, these 
matrix elements would not be all zero, and therefore the matrix elements Mj 4 (i — 1 — 3) 
would not be zero, which are responsible for the scalar-psuedoscalar mixing and hence CP 
violation. In short, these three CP phases in the radiative corrections generate the CP 
mixing between the scalar and pseudoscalar components, implying that the four neutral 
Higgs bosons are not states of definite CP parity. 
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The four physical neutral Higgs bosons at the one-loop level are defined as the eigen- 
states of the squared mass matrix M, with their squared masses as the eigenvalues of M. 
Let us denote the physical four neutral Higgs bosons as hi {i= 1 to 4) and their squared 
masses as m\, (i — 1 to 4). We assume that m\. < m\. for i < j. Among them, the upper 
bound on m\ x , the squared mass of the lightest neutral Higgs boson, may be obtained by 
noticing that the smallest eigenvalue of a positive symmetric matrix cannot exceed the 
smaller eigenvalue of its upper left 2x2 submatrix. Thus, the upper bound on m 2 ^ at 
the one-loop level is given as 

m 2 hl < \ 2 v 2 sm 2 2p + m 2 z cos 2 2p + 2g' 2 v 2 (Q 1 cos 2 p + Q 2 sm 2 p) 2 
3mt (Xscotf3A h -A t A i2 ) 2 2 2 
+ 8ttV (m 2 f -m 2 ) 2 9{ ^ ^ 



1 1 

mz m~ 



3m 4 b (Astan/3A Si - A b A h 



8ir 2 v 2 (m~ — m~ ) 2 



+ 



47r 2 w 2 (m? — m? , 

'fc 2 ' 



3m\\ 2 v 2 A| (yf(m 



87r 2 s 2 (m? — m~ ) 2 167r 2 s 2 




where the first three terms come from the tree-level Higgs potential while the other terms 
come from the one-loop corrections due to top quark, bottom quark, the exotic quark and 
their superparticles. 



III. NUMERICAL ANALYSIS 

For numerical analysis, we assume that the extra U(l)' emerges from the Eq group, 
as described in the previous section. To be concrete, we take the i/-model [19], where 
U(l)' is a mixture of U(l) x and U(l)^ with the mixing angle 8e = tan^ 1 y/lE. At the 
electroweak scale, the analysis of renormalization group equation leads to g\(vnz) — 0.46 
for U(l)' gauge coupling constant [19]. Then, the U(l)' charges of H 1 , H 2 , and S are 
given respectively as Q 1 w -0.4910123, Q 2 w -0.2995571, and Q 3 = -{Qi + Q 2 )- We set 
quark masses as m t = 175 GeV, m = 4 GeV, and rrik = 700 GeV. The renormalization 
scale is set as A = 700 GeV in the one-loop effective potential. We assume that the 
lighter squarks are larger than top quark mass. For the remaining parameters, we set the 
ranges for their variations as follows: 100 < A t = A = Ak,rfiA (GeV) < 2000, 100 < 
msusv = m T = m B = m K (GeV) < 1000, 1500 < s (GeV) < 2000, < <j) t AbAk < 2tt, 
1 < tan/3 < 30, and < A < 0.85. Note that one need not vary O since it appears 
always together with A\ in the definition of m a- Also note that the upper bound on A is 
chosen by considering that the the maximum value of A is about 0.83 in the analysis of 
the renormalization group equation. A general comment on s may be that a large s leads 
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to a heavy Z' in the present model. Moreover, the large s can satisfy the experimental 
constraint that the mixing angle between Z and Z' should be smaller than 2-3 xl0~ 3 . 

We calculate the upper bound on by varying all the relevant parameters within 
their allowed ranges for a given tan/3. The results are shown in Fig. 1, where the upper 
bound on is plotted as a function of tan /3, The solid curve is obtained by considering 
only the tree-level Higgs potential, whereas the dashed curve is obtained by considering 
the full one-loop Higgs potential with explicit CP violation. It is observed in the MSSM 
that the lightest neutral Higgs boson mass increases steadily as tan /3 increases. In the 
present model, where the MSSM is incorporated with an extra U(l)' , the lightest neutral 
Higgs boson mass does not increase in accordance with the increase of tan /3, because of 
both the Higgs singlet and Vd contributions due to an extra U(l)'. One can see that in 
Fig. 1 the lightest neutral Higgs boson mass becomes maximum at tan/3 ~ 1.5 for both 
curves. 

As we have remarked in the previous section, the scalar-pseudoscalar mixing, which 
is responsible for CP violation in this model, is triggered by three CP phases, <frt, 4>b, and 
<ftk- The size of the CP mixing between the scalar and pseudoscalar Higgs bosons could be 
regarded as maximal if sin0 4 = sin</>b = sin^ = 1, while there would be no CP mixing if 
sinc^f = sin^f, = sin^ = 0. In order to see a more clear picture about the amount of the 
CP mixing among the scalar and pseudoscalar Higgs bosons, we introduce (i = 1,2,3) 



Uitt*) = T77 , TT7 i > ( 19 ) 



2\M lA 




\M ti \ 


+ 


M 44 | 



which depend implicitly on other parameters as well. 

These Ui(4>k) measure by what amount <pk contributes to the size of CP mixing. In 
order to express the relative contributions by </>&, we further introduce 

Ui{<pk = 0) 

By evaluating f2j ((/>£. ), we may study the role of the exotic quark sector in the CP violation 
in the present model. 

In Fig. 2, we plot Q 3 as a function of <pk, where we fix tbsusy = — A t = s/2 = 1000 
GeV and A = 0.3, and we set for simplicity <p t = <Pb- The four curves in Fig. 2 correspond 
to four different sets of (f>t = <ftb and tan/3, namely, the solid curve is obtained for cf) t = 
06 = 7r/6 and tan/3 = 2, the dashed curve for <pt = 4>b = tt/6 and tan/3 = 20, the dotted 
curve for <p t = <p b = 7r/3 and tan/3 = 2, and the dash-dotted curve for <f> t = (f> b = n/3 and 
tan j3 = 20. It is trivial to notice that f2 3 = 1 at </>& = 0, n, or 2ir. For other values of </>&, 
we find that ^3 may change up to as much as about 5 x 10 3 for some parameter values. 
This implies that the effect of <pk upon Q3, or 003, is very significant. In other words, the 
matrix element IM34I depends critically on the complex phase of the exotic quark sector. 

We study the behaviors of Q± and ^2 for < </>fc < 2n, fixing the values of other 
parameters as in Fig. 2. Unlike f2 3 , they do not fluctuate widely against the variation of 
<pk. Thus, it can be said that the complex phase of the exotic quark sector has a strong 
influence on the magnitude of the mixing between /13 and h&, whereas it has relatively 
weak influence on the magnitudes of the mixing between h\ and or between \i2 and h^. 
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We calculate the masses of the neutral Higgs bosons for the same parameter setting 
as in Fig. 2: ttisusy = rriA = A t = s/2 = 1000 GeV, A = 0.3, together with four different 
sets of 4> t = <pb and tan/3, for < <pk < 27r. We obtain approximately 126, 1000, 1002, 
and 1043 GeV, respectively, for m^, rrih 2 , mh 3 , and m^. We find that these values are 
almost unchanged for < 0& < 2tt, and nearly the same within a few GeV for the four 
different sets of <p t = 4>t and tan j3. 

The mass of the extra gauge boson Z' is also estimated. We find that m^/ is stable at 
approximately 1039 GeV for the parameter values in Fig. 2. The mixing angle between Z 
and Z', denoted as olzz 1 -, depends on both tan/3 and s. We find that azz' = 1-39 x 10~ 3 
for tan/3 = 2, and and slightly increases for larger values of tan/3 as 2.93 x 10~ 3 for 
tan/3 = 20. 

We have studied the dependence of the mixing element between /i 3 and on the 
complex phase of the exotic quark sector. Meanwhile, one can evaluate effectively the 
absolute size of the CP violation by calculating the dimensionless parameter 



p^tAbAk) = ty\O u 21 31 41 \ , (21) 

where Oy are the elements of the orthogonal matrix which diagonalizes the squared mass 
matrix for the neutral Higgs bosons. The range of p goes from to 1 since the elements 
of the orthogonal matrix satisfy the orthogonality condition of J2j=i 0\ = 1- If P — 0, 
there would be no CP violation, whereas the CP symmetry would be maximally violated 
if p — 1. The maximal CP violation that leads to p = 1 takes place when 0\ x = 0\ x = 
2 31 = Ol = 1/4. 

In order to figure out the dependence of p on fc , we introduce 

PM = P^hM . (22) 
PVPt, <Pb, <Pk = 0) 

We calculated pk as a function of 4>ki f° r the same parameter setting as Fig. 2. If pk 
remains at 1, it would imply that pk does not depend on 0^. On the other hand, <pk would 
contribute more significantly to CP violation if pk moves farther away from 1. The result 
is shown in Fig. 3. The four curves correspond to the four sets of parameter values. Fig. 
3 shows that the fluctuation of pk is larger for tan /3 = 20 than for tan f3 = 2. This implies 
that (pk play an important role in p^ for large tan /3. One can see that pk fluctuates by up 
to 20 % for nonzero (pk- Since p(4>t, 4>b, 4>k) is an absolute measure for the CP violation in 
our model, it is reasonable to expect that the complex phase of the exotic quark sector 
may change up to 20 % of the CP violation. 



IV. SPONTANEOUS CP VIOLATION 

In this section, let us examine briefly whether it is possible to realize the spontaneous CP 
violation in the present model. At initial stage, the Lagrangian density of the model is 
assumed to be invariant with respect to CP property. After the electroweak symmetry 
breaking, the CP symmetry may spontaneously be broken by complex phases in the VEVs 
of the neutral Higgs fields. In the MSSM with an extra U(l)', one can conjecture from the 
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Georgi-Pais theorem that the spontaneous CP violation leads to a very light pseudoscalar 
Higgs boson, namely, the axion. It has been addressed that the spontaneous CP violation 
through radiative corrections can be realized when at the tree level there exist massless 
Higgs bosons other than Goldstone boson [23]. In the present model, there is no pseudo- 
Glodstone boson, which is referred to as axion, and thus the spontaneous CP violation 
cannot occur at the one-loop level. 

In order to examine the possibility of spontaneous CP violation scenario at the one- 
loop level, we can quantitatively analyze the axion mass by using the formulae which are 
derived for explicit CP violation scenario. In spontaneous CP violation scenario, the one- 
loop effective potential of the present model may have only one physical CP phase (0 S ) 
arising from the VEV of the Higgs singlet by considering the radiative corrections due 
to top quark, bottom quark, the exotic quark, and their superpartners. At the one-loop 
level, we obtain the minimum condition for the vacuum with respect to <p s as 

A\ = 3 li A \ f(m 2 f ,mj)+ l m ! A \ f(ml,ml) 
A 167r 2 t? 2 sin 2 /T V tiJ 167rVcos 2 /r V 6l ' h *> 

3m fc y4fc . 2 2 \ /no\ 

H 7rr:l{iTi7 , m r 23 

J v fci ' k 2 > \ ' 



for <p s 7^ 0, 7r . We obtain the axion mass at the one-loop level as 

2 3m 4 A 2 A 2 s 2 sin 2 4 g{m\, m|) 3m 4 A 2 A 2 s 2 sin 2 b 9^, mfj 



m 



8ir 2 v 2 sin 4 (3 cos 2 a (rn 2 2 — ) 2 87r 2 t> 2 cos 4 j3 cos 2 a (m~ — m- ) 2 
| 3mj\ 2 Alv 2 sin 2 fakrnlj ml) 
87r 2 s 2 cos 2 a (m~ — m|) 2 



(24) 



The massless function g in the above expression is found to have negative or very small 
values for the entire parameter space we consider, in the case of explicit CP violation 
scenario. Employing this result from the explicit CP violation scenario, the squared mass 
of axion is almost always negative for the entire parameter space, which is unacceptable. 
Therefore, the spontaneous CP violation is not viable for the MSSM with an extra U(l)'. 

Furthermore, the lightest neutral Higgs boson mass also faces difficulty in the spon- 
taneous CP violation scenario. Calculate the lightest neutral Higgs boson mass from 
the 4x4 squared mass matrix for neutral Higgs bosons, we find that the squared mass 
of the lightest Higgs boson has negative values for the above parameter space. In fact, 
the minimum condition for the vacuum with respect to cf) s has a very strong constraint, 
and thus the vacuum becomes unstable for the full parameter space we consider, in the 
spontaneous CP violation scenario. 



V. CONCLUSIONS 

We consider CP violation in the MSSM with an extra U(l)' symmetry, which is assumed 
to be originated from a string-inspired Eq model. The charge of the extra U(l)' is defined 
as Q' = cos 9eQ x + sin OeQ^p with 6e = tan -1 a/15. This model possesses one Higgs singlet 
and exotic quarks and squarks, in addition to the particle content of the MSSM. In this 
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model, the CP violation is impossible to occur either explicitly or spontaneously in its the 
tree-level Higgs sector. The mixing between the scalar and pseudoscalar Higgs bosons, 
which induce the CP violation, may occur via the radiative corrections at the one-loop 
level. For the radiative corrections, we consider the loop contributions due to top quark, 
bottom quark, the exotic quark and their superpartners. The radiative corrections contain 
in general three complex phases originated from the top squark sector, the bottom squark 
sector, and the exotic squark sector. We pay attention to the complex phase <pk that 
comes from the exotic squark sector. 

We calculate the lightest neutral Higgs boson mass. In explicit CP violation scenario, 
we find that the mass of the lightest neutral Higgs boson is smaller than about 135 GeV 
and 165 GeV, respectively, at the tree level and at the one-loop level, for the whole 
parameter space that we considered. 

For four different sets of parameters, we calculate the magnitude of the CP mixing 
among the scalar and pseudoscalar Higgs bosons for < <pk < 27r. The matrix element, 
I Mza | , representing the CP mixing between the heaviest scalar and the pseudoscalar Higgs 
bosons, is found to depend heavily on the complex phase from the exotic squark sector, 
4>k- In particular, the magnitude of the CP mixing between the heaviest scalar and the 
pseudoscalar Higgs bosons is found to change as much as about 20 % for < (fik < 27r 
when <pt = 4>b = tt/3 and tan/3 = 20, whereas the CP mixing among the other scalar 
Higgs bosons and the pseudoscalar Higgs boson is found to be relatively stable against 
the variation of 

Also, we have investigated the possibility of spontaneous CP violation at the one- 
loop level in the neutral Higgs sector, where the tree-level Higgs potential possesses CP 
symmetry. We argue that it is impossible to achieve spontaneous CP violation in the 
radiatively corrected Higgs potential of this model. 
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FIGURE CAPTION 



FIG. 1. : The upper bound on the lightest neutral Higgs boson mass as a function of 
tan /3, where other relevant parameters are allowed to vary within their respective ranges: 
100 < A t = A b = A k ,m A (GeV) < 2000, 100 < m SU SY = m T = m B = m K (GeV) < 1000, 
1500 < s (GeV) < 2000, < faAbAk < 2tt, 1 < tan/3 < 30, and < A < 0.85. The 
quark masses are set as m t = 175 GeV, = 4 GeV, = 700 GeV, and the renormal- 
ization scale is set as A = 700 GeV. The solid and dashed curves are respectively obtained 
from the tree-level and the one-loop level Higgs potential with explicit CP violation. 

FIG. 2. : Plot of Q3 as a function of </>&, for four different sets of <p t = <ftb an d tan/3: 
<t>t — 4>b — vr/6 and tan (5 = 2 (solid curve), <p t = 4>b = it / '6 and tan (3 = 20 (dashed curve), 
4>t — <fib = tt/S and tan/3 = 2 (dotted curve), and <p t = (f>b = vr/3 and tan/3 = 20 (dash- 
dotted curve). The remaining parameters are set as msusY = ttia = A t = s/2 = 1000 
GeV and A = 0.3. 

FIG. 3. : Plot of pk as a function of <pk, for the same parameter setting as Fig. 2. 
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FIG. 1: The upper bound on the lightest neutral Higgs boson mass as a function of 
tan f3, where other relevant parameters are allowed to vary within their respective ranges: 
100 < A t = A b = A k ,m A (GeV) < 2000, 100 < m SU SY = m T = m B = m K (GeV) < 1000, 
1500 < s (GeV) < 2000, < fa, <l>b, <l>k < 2tt, 1 < tan/3 < 30, and < A < 0.85. The 
quark masses are set as m t = 175 GeV, m& = 4 GeV, = 700 GeV, and the renormal- 
ization scale is set as A = 700 GeV. The solid and dashed curves are respectively obtained 
from the tree-level and the one-loop level Higgs potential with explicit CP violation. 



17 




FIG. 2: Plot of as a function of for four different sets of <pt = 4>b and tan/3: 
<t>t — <Pb = tt/6 and tan j3 = 2 (solid curve), <pt = <pb = tt/6 and tan (3 = 20 (dashed curve), 
0t = <Pb = tt/3 and tan/3 = 2 (dotted curve), and <pt = 4>b = tt/3 and tan/3 = 20 (dash- 
dotted curve). The remaining parameters are set as msusY = f^A = A t = s/2 = 1000 
GeV and A = 0.3. 
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FIG. 3: Plot of pk as a function of (pk, for the same parameter setting as Fig. 2. 
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